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Abstract. This paper investigates the properties of spin- 1/2 particle Zitterbewegung 
in the presence of a general curved space-time background described in terms of Fermi 
normal co-ordinates, where the spatial part is expressed using general curvilinear 
co-ordinates. Adopting the approach first introduced by Barut and Bracken for 
Zitterbewegung in the local rest frame of the particle, it is shown that non-trivial 
gravitational contributions to the relative position and momentum operators appear 
due to the coupling of Zitterbewegung frequency terms with the Ricci curvature tensor 
in the Fermi frame, indicating a formal violation of the weak equivalence principle. 
Explicit expressions for these contributions are shown for the case of quasi-circular 
orbital motion of a spin-1/2 particle in a Vaidya background. Formal expressions also 
appear for the time-derivative of the Pauli-Lubanski vector due to space-time curvature 
effects coupled to the Zitterbewegung frequency. As well, the choice of curvilinear co- 
ordinates results in non-inertial contributions in the time evolution of the canonical 
momentum for the spin-1/2 particle, where Zitterbewegung effects lead to stability 
considerations for its propagation, based on the Floquet theory of differential equations. 



PACS numbers: 03.65.Ta, 04.20.-q, 04.20.Cv 
1. Introduction 



The concept of Zitterbewegung, or "trembling motion" in German due to the rapid 
oscillation of a spin-1/2 particle about its classical worldline, is one that still generates 
great interest about foundational issues in quantum mechanics, even after over 70 years 
since its initial discovery by Schrodinger [1]. When investigating the free-particle motion 
of an electron via the Dirac equation, he showed that its instantaneous velocity winst. 
equals the speed of light c, which clearly disagrees with the fact that a massive particle's 
observed speed will always be fobs. < c. However, Dirac convincingly emphasized [2] that 
this apparent contradiction is resolved by noting that fobs, is really an average velocity 
determined by measuring the electron's position over a sufficiently small time interval 
At. When cAt ^ h/mc, the Compton wavelength for the electron of mass m, it follows 
that the practical measurement of fobs, always leads to a preservation of causality, and 
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without any complications due to Zitterbewegung. Nevertheless, this highly non-intuitive 
oscillatory behaviour identified to exist for massive spin-1/2 particles is a conceptually 
rich feature of rclativistic quantum mechanics that still offers both thought-provoking 
challenges and opportunities for deeper exploration [3, 4, 5], especially when applied to 
more complicated situations. 

While past and recent investigations of Zitterbewegung since Schrodinger's and 
Dirac's studies have been largely motivated by attempts to understand the intrinsic 
structure of the electron [6], there currently exists an intensive practical study of 
the concept within the context of condensed-matter theory and experiment. Recent 
explorations of this phenomenon involving two-dimensional carbon sheets known as 
graphene have been proposed [7, 8], where Zitterbewegung effects are potentially 
observable in the non-relativistic limit of the Dirac equation, while a second paper 
claims a direct observation for photons in a two-dimensional photonic crystal [9], and 
a third paper proposes detection of Zitterbewegung using ultracold neutral atoms [10]. 
These studies are especially interesting given the largely accepted understanding that 
free-particle Zitterbewegung arises only in situations where positive- and negative-energy 
solutions to the Dirac equation can interact, allowing for quantum interference effects 
to manifest the rapid oscillation terms that would otherwise remain decoupled and 
unobservable. 

This leads to the following consideration where it concerns interactions in a 
gravitational background. Given that external fields can introduce non-trivial effects 
that lead to observable consequences for electron Zitterbewegung, and given that 
gravitation is best described (for now) by Einstein's theory of general relativity in 
terms of the dynamical curvature of space-time due to external matter, it is reasonable 
to envision curved space-time as something like a medium with an energy density 
content not unlike that of an electromagnetic background. Of course, there is the 
fundamental distinction that Einstein's equivalence principle is an integral part of 
general relativity, rcficctcd mathematically as diffeomorphism invariance under co- 
ordinate transformations. This is certainly reasonable when considering purely classical 
descriptions of matter, where it is meaningful to speak about its localization within a 
compact spatial region. However, conceptual difficulties arise when considering quantum 
matter in a curved space-time background, where issues of non-locality due to the wave- 
particle duality of quantum objects overlap with classical descriptions of space-time 
curvature. This is a particularly relevant issue considering that gravitational effects are 
manifested locally from observing tidal forces via the geodesic deviation equation, which 
requires the existence of a neighbouring geodesic in relation to the reference worldline. 
It is very unclear how to disentangle the competing effects of quantum non-locality and 
the locality of classical gravitation in a way that results in meaningful statements about 
quantum matter effects in a curved space-time background, including Zitterbewegung. 

The purpose of this paper is to explore the physical properties of spin-1/2 
particle Zitterbewegung in a curved space-time background described locally in terms 
of Fermi normal co-ordinates. Its motivation is based on understanding whether spatial 
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fluctuations about the particle's worldline lead to non-trivial gravitational corrections 
that influence its propagation in space-time. A recent innovation in the standard 
formalism of Fermi normal co-ordinates is to describe the spatial fluctuations in terms 
of curvilinear spatial co-ordinates [11, 12, 13]. This paper begins with a brief review 
of the existing treatment of Zitterbewegung in a flat space-time background, as given in 
Section 2. Because the Fermi normal co-ordinate system is defined locally about the 
particle's worldline, it is important to make use of a local perspective for the occurrence 
of Zitterbewegung and explore its physical consequences. It so happens that the approach 
introduced by Barut and Bracken [14] is particularly well-suited for this purpose, and so 
this perspective is adopted for consideration in a curved space-time background. This 
leads to a brief introduction of the covariant Dirac equation in Fermi normal co-ordinates 
in Section 3, which includes the curvilinear co-ordinate generalization. A presentation 
of the corresponding Dirac Hamiltonian derived from within this framework is given 
in Section 4, where contributions due to Zitterbewegung reveal a formal violation of 
the weak equivalence principle due to direct coupling of the gravitational field with 
the particle's mass. From this Hamiltonian, it becomes possible to compute the time- 
evolution of quantum operators in Section 5 and show how Zitterbewegung formally 
appears in the expressions. Particular attention is given to the time-evolution of the 
position and momentum operators, and also the covariant spin operator described by 
the Pauli-Lubanski vector. A brief conclusion follows in Section 6. 

For this paper, geometric units of G = c = 1 are adopted, while the curvature 
tensor definitions follow the conventions given by Misner, Thorne, and Wheeler [15], 
but with metric signature —2. As well, the flat space-time gamma matrices follow the 
conventions of Itzykson and Zuber [16]. 



2. Review of Free-Peirticle Spin-1/2 Zitterbewegung in Flat Space-Time 

2.1. General Formalism 

Before introducing Zitterbewegung in curved space-time, it is useful to first review the 
original treatment of the problem in flat space-time by Schrodinger [1, 2, 14]. Assuming 
space-time co-ordinates in Cartesian form X^^ — {T,X), the Schrodinger equation is 
ih (d/dT) iP{X) = HoiP{X), where 



Ho = m/3(o) + a(o) 



(0) 



(1) 



is the free-particle Dirac Hamiltonian defined in terms of spin-1/2 particle mass m, 
canonical momentum P^°^ = ihd/dX^, and 4x4 matrices a.L^ — 7° 7-' and /3(o) = 7° 



with I7" ,7'^! = 277"^ [16]. As expected, the canonical momentum operators satisfy 



rp(o) p(c 



(0) 



0, 



while the corresponding position operators XL. satisfy 



^ (0) ) (0) 



0, 



(0) 

p(o) 

^ (0) > 



-ihS' 



(2) 



(3) 
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Adopting the Heisenberg picture, it follows that the time derivative of some operator 



A is 



dA(r) 



[^0 , A{T)] 



(4) 



dT h 

Clearly, dP-'^^-^{T)/dT = dHo/dT = 0, while the instantaneous velocity and acceleration 
operators are 



dT 



- -^'n^Q^[i ) Mo, 



dT^ 



dT 



where 



'(0) 



(0)^ 



(5) 
(6) 
(7) 



Because P|q^ (T) and Hq are constants of the motion, it is also true that 



dr7;o)(^) 
dT 



dT 



PUT) Ho \ 



dayr) 

dT 



2i 



(8) 



- -^^io)(T)Ho. 
An immediate solution follows for r7|g^ (T) , such that 

or 

c^UiT) - PUiO) H^' + r7j;)(0) e-^^W^ . 

Direct integration of (10) then leads to the free-particle solution for X|q-j(T), such that 
the time-evolution for the position and canonical momentum operators are 

x\o){t) = xUq) + {pUq)h^')t 

(11) 



(9) 
(10) 



+ I r7jo)(0) [sin {2HoT/n) - 2i sin^ (HoT/h) 



PU(T) =P?o)(0). (12) 

The first two terms of (11) denote the expected mean position as a function of time T, 
while the remaining terms show the rapid oscillation effect superimposed about the spin- 
1/2 particle's worldline due to Zitterbewegung. For future reference, it is straightforward 
to show that Zitterbewegung also applies to P(o){T) and the Pauli spin matrix <t|q)(T), 
such that 



Ao)(T) 



= (3o e-^'"°'^/^ + m [2 sin^ (HoT/h) + i sin (2HoT/h)] , 



(13) 



a^T) ^(rl + i (ao x P(o)(0))' H^' [2 sin^ (HoT/h) + i sin {2HoT/h)] , (14) 
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where 

a^o = aj;)(0), /3o = /3(o)(0), = ctI,^{0) . (15) 

Before continuing, it is worthwhile to explore the mathematical details of the 
operators which describe the Zitterbewegung phenomenon. Based on a precisely 
formulated description of the problem [17], it is possible to establish that the domain T> 
corresponding to the position operator X|q^ (0) is invariant under time evolution, such 
that 

V (Xfo) (T)) = e-'"°^/^ V {X{,^ (0)) =V (Xfo) (0)) . (16) 

As such, (16) implies that X-^q^(T) remains bounded for all finite T. In addition, 
the expectation value for the time-evolved position operator is well behaved, in that 
the component of (11) exhibiting Zitterbewegung integrates to zero in the expectation 
value as T — > oo. Furthermore, it is well-known [17, 18, 19] that the Hilbert space 
identified with solutions of the free-particle Dirac equation decouples into invariant 
subspaces defined by positive and negative energy states, respectively. This implies 
that the expectation values of ^(o) {T) formed by exclusively positive energy or negative 
energy states never reveal any Zitterbewegung-induced effects. Not only does it explain 
why Zitterbewegung is not an observable for a propagating free particle, it suggests the 
existence of a breakdown in the single-particle treatment of quantum matter, justifying 
the need for the standard quantum field theory description. 



2.2. The Barut-Bracken Approach to Zitterbewegung 

The main expressions above can apply to spin-1/2 particle motion in a general reference 
frame with P|q^(T) 7^ 0. Schrodinger's original view of Zitterbewegung was to search for 
"microscopic" degrees of freedom for position, momentum, and angular momentum to 
explain the electron's rapid oscillatory behaviour as a separable effect superimposed on 
its observable motion [1, 14]. To this end, he proposed the existence of a microscopic 
position operator ^-'(T), according to 

Xl,){T)^Xl^{T)+e{T), (17) 



where 



X{a) (T) = Xfo) (0) + (p|o) (0) H^') T - y ryfo) (0) H^' , (18) 



ih 

along with a microscopic momentum associated with ^■'(T), in the form 



m =v<(0)^o-^e-^^W^, (19) 



dXUT) 

rjl,^ (T) H, + (T) = Ho . (20) 

In re-examining Zitterbewegung with the motive to understand the electron's 
intrinsic structure for addressing sclf-cncrgy and renormalization issues in QED, Barut 
and Bracken [14] saw no difficulty with Schrodinger's formulation of microscopic 
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position, according to (17). However, they had serious issues with his definition 
of microscopic momentum (20), primarily motivated by the claim that an equally 
attractive-but incompatible-definition could be proposed within his formalism, such 
that 

Ho r){,) iT) + Pfo) (T) = H, . (21) 

Another major deficiency is that (20) and (21) are not Hermitian operators, and that 
it is impossible to construct a viable Hermitian combination of these two momentum 
definitions. As well, because the canonical momentum P|q-)(T) is a constant of the 
motion, they argue that it should represent the total overall momentum of the system, 
and that the microscopic nature of (20) should not have any dependence on an inherently 
macroscopic description of momentum. 

To deal with these conceptual deficiencies identified by Barut and Bracken in 
Schrodinger's treatment of Zitterbewegung, they proposed the existence of a relative 
momentum operator 7'(g-)(T) that results from performing the computations in the 
local rest frame. || This amounts to letting the operators A(o)(T) — > A(b)(T), where 

P{b){T)^Q. (22) 

By also setting the initial position to X|g^(0) = 0, a relative position operator in the 
local frame can be identified, such that 

X\g^ (T) ^ ^ [sin (cu.itt.T) + 2i sin^ {uj^,,T/2) /3o] a'o , (23) 
while the remaining operators (10), (13), and (14) become 

ajg) (T) ^ al [cos (c^zitt.^) - i sin (u;.itt.r) po] , (24) 

PiB){T)^(3o, (25) 

al^^{T)^al (26) 

where cUznt. — '^'m/h is the Zitterbewegung frequency derived from 

Ho =m/5(B)(T) ^ m/3o: (27) 
H^'^m-'f^o. (28) 

Barut and Bracken subsequently propose that the relative momentum is described by 

Pj^)(r)=maJ^)(r) (29) 

with respect to the spin-1/2 particle's classical worldhne. 

It should be noted that X^(^b) (^) differs slightly from ^^{T) incorporated in the local 
rest frame by Barut and Bracken [14]. This is due to a minor change in definition, such 

II The main advantage of Barut's and Bracken's approach over Schrodinger's is that the derived relative 

position and momentum operators form a harmonic oscillator relationship with frequency Wzitt. and are 
part of a basis set satisfying an S0(5) Lie algebra structure [14]. 
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that the third term of (18) is absorbed by (19) to form (23), instead of expressing (11) 
as given by (17)-(19) in the local rest frame, where 

iti 
2m ' 



(30) 



ih 
2m 



h 

2m 



[sin (cUzittT) - i cos (cUzutT) Po] 



^ J'^.itt.(T--K/2)l3o 3 



(31) 



With the definitions for relative position and momentum now given, it is possible 
to derive the uncertainty in their measurements, to compare with the Heisenberg 
uncertainty relation. To begin, consider the time-averaged operator over a complete 
cycle in terms of LO^y^x, , such that 



O^zitt. 



A{T) dT . 



27r h 

It follows that the Zitterbewegung time average of (23) is 

ih 



and subsequently leads to 



1/2 



73 /^^ 
2 \m^ 



(32) 
(33) 
(34) 



in agreement with an earlier computation derived in a similar context [20]. However, 
the squared- magnitude of (T) prior to time averaging is 



12 



sin^ (a;zitt.r/2) , (35) 



Zltt. 



which leads to 



1/2 



h 



, m 



Therefore, the local uncertainty in relative position is 
AX 



■(B) 



X<B){T) 



X<B){T) 



2 [m, 



(36) 



(37) 



To compute the corresponding local uncertainty in the relative momentum, it is 
self-evident that 



while 



rU{T)) =0, 

2\ 1/2 

^(B)(T) 



m . 



This results in ^'P(b) — V^m, leading to the expression 

well in agreement with the Heisenberg position-momentum uncertainty relation 
(AX(b)) (AP(b)) > h/2. 



(38) 

(39) 

(40) 
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3. An Approach to the Curved Space-Time Generalization 

3.1. The Hypothesis of Locality 

Fundamental to the requirements of this paper is consideration of the imphcitly accepted 
hypothesis of locality [21] for classical phenomena and how it relates to the intersection of 
quantum mechanical behaviour in curved space-time. This hypothesis essentially states 
that, at any instantaneous moment of proper time on a classical worldline, it is possible 
to identify an instantaneous comoving ineriial frame tangent to the worldline. A smooth 
one-to-one correspondence follows from this identification, along with an intrinsic length 
and time scale correlated with the comoving observer. 

When considering purely classical matter, this hypothesis is undeniably successful. 
Issues become much more ambiguous, however, when applying this hypothesis to 
single-particle quantum states-let alone quantum fields, since the wave-particle duality 
implies the need for a region of space-time to establish the location of quantum matter 
that is consistent with Heisenberg's uncertainty principle. This especially applies to 
quantum matter with long wavelengths when compared to the characteristic length of 
the background gravitational source, where quantum interference effects become most 
relevant to disentangle conceptually. If a measurement applied to quantum matter 
is somehow performed on a sufficiently small time scale compared to the comoving 
observer's intrinsic time scale, it is reasonable to surmise that the hypothesis of locality 
still applies. However, this claim is still a tenuous one, at best, and should be subjected 
to deeper analysis. 

3.2. Covariant Dirac Equation in Fermi Normal Co-ordinates 

Accepting the hypothesis of locality as a viable starting point, the approach to 
Zitterbewegung in curved space-time is to adopt Fermi normal co-ordinates with 
respect to a local comoving frame [11], with X° = T identified as the spin-1/2 
particle's proper time and as the Cartesian spatial co-ordinates orthogonal to 
the worldline. Again, for purely classical phenomena, this choice for X^ is perfectly 
reasonable. However, when considering the propagation of quantum matter, it may be 
more appropriate to use general curvihnear co-ordinates U*^ — (T, u^) to better reflect 
the symmetries associated with the particle's trajectory in space for general motion 
between neighbouring intervals of proper time, where X^ — X^{u) [11]. Eventually, the 
spatial Fermi normal co-ordinates will be identified with quantum fluctuations about 
the classical worldline, such that X^ X|^-j(T) when deriving the Dirac Hamiltonian. 
This is the central hypothesis that results in the description of curved space-time 
Zitterbewegung to follow. 

Starting with a worldline C defined in a general space-time background and 
parametrized by proper time r, the Fermi frame [22, 23, 24, 25] is determined at some 
event Pq on C by constructing a local orthonormal vierbein set {A"^} and inverse 
set {X^a}- Then the local spatial axes are defined by A'^o — dxP'/dr and A'^^. If 
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^P- — {dxP'/da)^^ denotes the unit spatial tangent vector from Pq to a neighbouring 
event P, where a unique spaceUke geodesic orthogonal to C exists along proper length 
a, then T = r and X* = a^^X^p, become the Fermi normal co-ordinates at P. The 
corresponding space-time metric is then described by 

ds^^^g,,{X)dX^dX\ (41) 



where 



^^oo(X) = 1 + ^Roi^AT) X'X^ + --- , (42a) 
^^oi(^) = \ ^Ro^Jk{T) X^X' + ---, (426) 

^gij{X) = r]ij + ^ ^RikjiiT) X" X' + ■ ■ ■ , (42c) 



and 



e°,•(X) = l^i?o,,■fe(^)X^X^ (45&) 



^Ral3-fS{T) — Rflupa X^a X" i3 X'^ S ■ (43) 

The covariant Dirac equation for a spin-1/2 particle with mass m can be written 
in terms of (42a)-(42c) as 

[ir{x) {d, + i r,{x)) - m/h] i;{X) = , (44) 

where 9^ = d/dX^^ and r^(X) is the spin connection defined with respect to {7''(X)}, 
the set of gamma matrices satisfying {'y'^{X),Y{X)} — 2 gF^i^)- ^ local Lorentz frame 
[26] can be determined according to ^g'^i,(X) = 77^^ e"^(X) e^i,{X), where {e'^a{X)} and 

|e°^(X)| form a respective orthonormal vierbein and inverse vierbein set, such that 

e%{X) = l + \ "^RoioAT) X'X^ , (45a) 
1 
6 

e'o{X) = - ^ '^R'MT) X'^ , (45c) 
e'jiX) = 5' J - ^ ^R\jiiT) X' X' , {A5d) 

e%{X) = 1 - ^ ^Roioj{T) X'X^ , (46a) 

e%{X) =^^R'MT)X^X\ (466) 

e°,-(X) = - ^ ""RoijkiT) X' X^ , (46c) 

eVX) = 5', + ^ ""R'kjiiT) X*^ X' . (46d) 
The spin connection is then determined to be 

= I nX) (V^7a(^)) - e\ (V^ e-^«) , (47) 



and 
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To first-order in 



where V^^ is the covariant derivative operator and (7°^ = | 
the Riemann tensor, it is shown from (47) that 



7°, 7^ 



ro(X)=2 7°7'' 



2 ^ RmkiT) + - ^ RijQk,o{T) 



X^ , 



Ti{X) =i7°y 

where ^Ri[ko]i{T) 
two indices. 



Rmk{T) + ^Ri\kQ\i{T) — — ^Rijki^o{T) X^ 



X\ 



(48) 
(49) 



Rim{T) — ^ RiQki{T) denotes antisymmetrization of the middle 



3.3. Conversion to Curvilinear Co-ordinates 

It is straightforward to introduce a conversion of Fermi normal co-ordinates from locally 
Cartesian to general curvilinear co-ordinates, in terms of a new set of orthonormal 
vierbeins described by 



e^a{U)^^e-a{X), 



dip 



e"a(X) 



(50a) 
(506) 



The corresponding spin connection in curvilinear co-ordinates becomes 
r.(C/) = J r{U) [V,7a(C/)] = (Vm e\) 



4 

dX'^ 
dip 



TaiX) . 



(51) 



It follows that the covariant Dirac equation (44), expressed in curvilinear co-ordinates 
and projected onto the local Lorentz frame, is 

V (Va + i r^{U)) - m/h] iP{U) = , (52) 

where = -\- itf^{U) is the flat space-time covariant derivative operator in 
curvilinear co-ordinate form. The corresponding line element is represented by 

ds^ = dT^ + rivj (a(*~) {u) du') (\^^^ {u) du>) , (53) 

where X^^\u) are dimensional scale functions [27], and 

d _ Id 



(54) 



A straightforward computation from (52) leads to the covariant Dirac equation in 
the form 



where 



r{P^c-hT^{U))-m\i;{U)^0, 



(55) 

(56) 
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is the canonical momentum operator in curvilinear co-ordinates, with 



(S) 

6 

(S) 

i 

(T) 



ih 



V^ln [\^''\u) \^^\u) X^^\u)) 



1/2- 



M ^ / / Of/ml " 



1 



0[/m] 



-'^R"'jmkfiiT)X^X', 

2 RjOOm{T) + - RjlOm,o{T) X 

Yri^,u{T)XK 



X' 
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(57) 

(58) 
(59) 
(60) 
(61) 
(62) 



For example, it is straightforward to verify in spherical co-ordinates — {r,9, (p) that 
[12] 



\dr r I ^ r \d9 2 /' 



p<t> — 



ih d 



r smi 



(63) 



where \'^^\u) = 1 , X^'^\u) = r , and X^^\u) = r sin 9 . 
The final step follows by recalling the identity [28] 

Y Y Y = rt'" 7^ - 2 -IS^'^ -iYY ^^^''a , (64) 
where ef^^P^ is the Levi-Civita symbol with e^^^s _ Use of (64) leads to a new 

expression for the spin connection, such that 

Ta =l'^f-^^t\ (65) 

=^''%fgL]> (66) 
where the "C" in (66) is the chiral-dependent part of the spin connection, while the "S" 
in (60)-(61) denotes the symmetric part under chiral symmetry. 



4. Spin-1/2 Particle Dirac Hamiltonian for General Motion in Curved 
Space-Time and Zitterbewegung Effects 

Having obtained the required computations for the covariant Dirac equation in Fermi 
normal co-ordinates, it is straightforward to determine from the Schrodinger equation 
ih (d/dT) ijj(U) = Hi/jiU) that 



H = ml3 + {a-P)+h 

ih 



a^B,{T,X)--a^£,{T,X) 
ih 



(67) 



is the Dirac Hamiltonian expressed in terms of the "electric'' and '"magnetic'' field 
components [29] 

S,{T,X) = -''R^^^k{T)X\ (68) 
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Bj(T, X) = ^ • ^RimokiT) X" , (69) 

where dj B^(T, X) = and dj {T, X) = —^Rqq(T), the analog of Poisson's equation 
in electromagnetism, whose matter field source generates the Ricci curvature tensor 
in the Fermi frame. The Dirac Hamiltonian (67) is comprised of operators which are 
defined in the local Lorentz frame and are related to corresponding operators in general 
curvilinear co-ordinates in the standard way via the vierbein projection functions (50a) 
and (50&). As such, the mathematical properties of these operators locally correspond 
precisely to the description given for their counterparts in strictly fiat space-time, with 
the understanding that the range of curvature deviation satisfies the condition that 
^ g^y{X) — rj^i, <^ 1, particularly when — » X|^^(T). Apart from the spin-dependent 
term coupled to (69), it is evident that the gravitational contributions to (67) in its 
present form are anti-Hermitian. If the time variations of the gravitational field are 
deemed small, then the last two terms in (67) can be neglected. However, for this paper 
all terms in the Dirac Hamiltonian are retained. 

At this point, it is important to note that the Dirac Hamiltonian H — > Hq — 
m(5 -\- {oL ■ P) as X^ — > 0, satisfying the weak equivalence principle for purely classical 
conceptions of space-time. However, it so happens that if X^ JC|g^(T) in the 
Fermi frame, then a quantum violation of the weak equivalence principle emerges due 
to Zitterbewegung. This can be shown explicitly by the following process, where for the 
sake of notational convenience the local Lorentz frame indices will be unhatted for all 
subsequent computations. Suppose the Dirac Hamiltonian is expressed as 

H = Ho + Hg, (70) 

where Hq is given by (1) and Hq is the interaction Hamiltonian involving space-time 
curvature. Given that \X^\ ~ h/m and the curvature tensors in the Fermi frame are 
defined on the classical worldline to satisfy ^gn^{X) — rj^^ <^ 1, it is reasonable to 
surmise that Hq <^ Hq. If X^ — > X|g^(T) in the Fermi frame as operators, then the 
X^ X"^ within (67) become 

X[b) (T) X^b) (T) = 1 {x;^) (T) , X^^) {T)} + 1 [X[s^ (T) , X^^) (T)] , (71) 

the sum of anticommutator and commutator expressions, respectively, involving 
X\s^(T) and X^)(r). This leads to 

= ^e''-/Ri^o,{T)crl^^{T)Xf^^{T) - | ^i?o,o.(T) aJ^)(T) Xf^)(T) 

- f ^RjiimMT) c^^B)iT) {x\B)iT) , X^^)(T)} 

+ ^ ""RjoimAT) Ix\b){t) , x^^)(r)] 

- I ^R'l^n^AT) {X[^){T) , X^B)(r)} (72) 

from substitution of (71) into the interaction Hamiltonian, with ^Rj(^im)o{T) = 
RjimoiT) + ^ Rjrnio{T) denoting symmetrization of the middle two indices, and 
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using the cyclic permutation property of indices for the Riemann tensor in the fourth 
term of (72). Upon evaluation in terms of 



{x;^)(T),X^^)(T)} = - -I- sin^ (a;,Ht.T/2) r^'™ 



zitt. 



8i 



sin' (u;zitt.T/2) e 



Olm 



n "0 ) 



(73) 
(74) 



zitt. 



and with substitutions of (23), (24), and (26) into (72), it follows that 

n 



G 



'^zitt. 



+ 2 



sin' (u;.itt.T/2) {^i?oo(T) 
4 



+ 



+ 



ih 



'^o,(r) + 

•J^zitt. 

sin (cj^itt.T) ^ 



(^zitt 

sin [uj^itt 



sin {u^ittT) ^Roj^o{T) 



3 l^zitt. 



sm u. 



zitt 



T/2) ^W',,,o{T) 



sin' (a;.itt.T/2) cos {u^^t.T) ^Rojfl{T) 



< , (75) 



3 cc;zitt. 

where the quantum violation of the weak equivalence principle formally appears due 
to the explicit coupling of space-time curvature to the spin-1/2 particle's rest mass via 
the Zitterbewegung frequency uj^itt. = 2m/h. This is by no means the first time that a 
suggestion of weak equivalence violation has appeared in the literature. For example, 
it was suggested years ago by Greenberger that, for a bound-state quantum mechanical 
particle in an external gravitational potential, various violations of the weak equivalence 
principle emerge due to the direct coupling of particle mass to the gravitational field 
that subsequently do not appear in the classical limit [30]. 

Clearly, (75) vanishes as /i — or m — oo, while the Hermitian part of Hq — >■ 
and the anti-Hermitian part of Hq is regular as m — > 0. In fact, it is shown that 

1 

12 

^i?o,(T) + ^^i?o,,o(T)T 



lim Hq 



ih 



T. 



(76) 



which is relevant for (almost) massless neutrinos. It is particularly interesting to note 
that the first two terms in (75) coupled to /5o form a gravitationally-induced efi^ective 
mass, where ^Rqq{T) is likened to a potential energy term, while the remaining terms 



|«o| 



1. A 



serve as energy flux contributions with instantaneous velocity Viast. 
similar identification exists for the anti-Hermitian terms not coupled to /3o, which yield 
decay width contributions due to Zitterbewegung in curved space-time. 

It is worthwhile to make a general comment about the anti-Hermitian contributions 
to (75). Since it is practically true that the time scale associated with the Fermi frame 
curvature tensor will be very much longer than the time scale corresponding to the 
Zitterbewegung frequency, it is possible to regard ^R^iuapiT) as effectively constant in 
time over a complete cycle defined by o^zitt.- When this is taken into account, (75) 
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reduces to 



Hg^ sin^ Kitt.T/2) {^i?oo(T) 



+ 2 



-^0^(7") + ^ sin {uj^itt. 



+ sin^ (^zittT/2) ^i?^^fe,o(T) 

- 8 cos (a;,ittT) ^i?o,-,o(T) a^o] ■ (77) 

If it so happens tiiat tiie curvature tensors time-evolve adiabatically [31] and ocqI3q — >• 
I |q:;^ , /9o| = 0, then the sole contribution to Hq is Hermitian. Therefore, in a practical 
sense the anti-Hermitian terms in (75) are of little consequence in determining the 
dynamics of a massive spin-1/2 particle in curved space-time due to Zitterbewegung. 

5. Time Evolution of Quantum Operators 

Having now determined the Dirac Hamiltonian (70) with Zitterbewegung contributions 
according to (1) and (75), it is possible to compute the time evolution of quantum 
operators via the Heisenberg equations of motion. This can be done in perturbative 
form by solving for A A{q) -\- , where 

dA(o)(T) ^ I 



dT n 



[^o,A(r)], (78) 



^-^^{[Ho^MT)]. (79) 

By this approach, explicit computations for the time evolution of the relative 
position and momentum operators can be determined in a gravitational background, 
where the specific case of quasi-circular motion around a spherical black hole described 
by the Vaidya metric is considered in detail [32]. This eventually leads to additional 
contributions to the position-momentum measurement uncertainties obtained in (40) 
due to space-time curvature and Zitterbewegung. A separate computation for the 
time evolution of the canonical momentum P-'{T) P|q-j(T) + P-'f^-y^{T) can also be 
determined, noting that in curvilinear co-ordinates the zcroth-ordcr terms P|q^(T) are 
not constants of the motion. When coupled with the periodic nature of the resulting 
vector differential equation, it follows that solutions to P\T) take the form given by 
Floquet's theorem [33, 34]. Finally, a formal presentation of the time derivative for the 
Pauli-Lubanski vector, the operator that describes particle spin in covariant form, is 
given with respect to the perturbation Hamiltonian Hq and briefly discussed. 
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5.1. Time Evolution of the Relative Position and Momentum in Curved Space-Time 
due to Zitterbewegung 

5.1.1. Relative Position X^{T) In the local rest frame, the relative position operator 
is 



where X|^-j(T) is given by (23) and 



dT 



CO. 



sin^ (a;,itt.T/2) 



zitt. 



X r^oo(T) 



sin^ + ^ sin {uj^itt.T) (io-i (5^ 



— I 



- 2 



'i?°^o(r) sin^ (a;,itt.r/2)/3o 



'Rokfi{T) sin (cUzitt.T) 



3 uj^itt. 

Retaining only the Hermitian part of (81) and integrating then leads to 

"^zitt. "'O 



'Roo{T') sin2 (a;,itt.T72) 



X 



Diagonahzation of (82) then results in 

4 . rT 



diag. 



, ,2 "0 



\ r^RooiT') sin\uj^sxt.T'/2)dT' , 
Jo 



where the squared magnitude of X^-, (T) before time averaging yields 



Xii){T) 



- r),j X\,^iT) Xl^{T) 

4^3 rT 



f ^i?oo(r')sin3(^,i„r/2)dT' 

<^zitt. "^0 



(80) 



(81) 



sin^ (a;.itt.T72) cxl + - sin {u^^t.T') (ioci Po) dT' . (82) 



(83) 



(84) 



5.1.2. Relative Momentum T'-' {T) With the relative momentum described by 

VKT)^moc^{T) ^ Vl^)iT) + Vl^{T), 
the time derivative of the first-order perturbation is 
dPji)(T) _ 2m 



dT 



+ 



'^zitt. 

2im 

'i^zitt. 

2im 

<^zitt. 



J^{T) [sin (u^itt.T) ai + i cos (u^iu.T) (3o 
cos (cu,itt.T) C^ (T) - sin (cu^HtT) X»^ (T)] /3o 
[sin (a;,itt.r) C\T) + cos (a;,itt.r) X>*=(T) 



'0 > 



(85) 



(86) 
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(87) 



= - 2 sin" (cu,itt.T/2) 



'RoAT) + 



+ 



3 uj^itt 

s\v? (a;zitt.T/2) cos (c^zitt.T) ^Rojfi{T) 



sin {uj^itt.T) ^RojfliT) 



sin {uj,^,T)^RoAT) 
8 



dT', 



3 cUzitt. 

By also retaining only the Hermitian part of (86) and integrating, it is shown that 

where in diagonalized form (90) becomes 

2m ■ rT 

riAT) 



1^ J^{T') [sin (a;,itt.r') a^o + cos (c^zittT') (i«'o A 



diag. cJzitt. " -'O 



'Rm{T') sin^ (u;zitt.r72) dT' . 



^9) 



(90) 



(91) 



A corresponding computation for the squared magnitude of 'P(i)(T) prior to time 
averaging then results in 



2^3 



m 



'i?oo(T') sin^ (a;,Ht.T72) dT' 



(92) 



5.1.3. Position- Momentum Measurement Uncertainties in the Vaidya Background To 
illustrate the contribution of Zitterbewegung in curved space-time to the position- 
momentum uncertainties in measurement, consider the quasi-circular motion of a spin- 
1/2 particle around a spherically symmetric black hole in the presence of radiation, 
as described by the Vaidya metric [32]. For this computation, it is assumed that the 
background source M monotonically changes in time according to 

M = Mo + AM(T) , (93) 

where Mq is the (static) Schwarzschild mass, 

d (AM) 



a 



AM(T) = - 



A" 



= 1 



2Mo 



T, 



(94) 
(95) 



r is the particle's orbital radius, and |d (AM) /d^| is the source's rate of change along ^, a 
radial null co-ordinate. If a = -|-1, then ^ is the advanced null co-ordinate corresponding 
to infalling radiation, while a = —1 implies that C, is the retarded null co-ordinate for 
outgoing radiation. Then it is determined that [32] 



"RooiT) 



2a 



d(AM) 



1 + 2a 



d(AM) 



dC 



X 



^"^^^.C(r,QKr)-^sin^(QKT) 



2 Are 



(96) 
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3Mo 



(r > 3Mo) 
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(97) 



and 



C{r,QKT) = 2 sin(2QK7') 
N 



+ — — [(1 - 2 tQk) sin (2 Qk^) - 2 Qk?'] 
To leading order in |d (AM) /d,^|, it is straightforward to show that 



16 



a 



(1)^ Vvaidya ~ 3 U^^,,^ 



d(AM) 



\X{T)\' 



Vaidya 



+ 



16 a 



zitt. 



d(AM) 



(98) 



(99) 



- — 7^ + i/?o7o 

OTT 



l(^(^))lvaidya ^ 

which results in 

(^-^)vaidya 



+ 



16 a 



d(AM) 



a 



256 



157r iV2r2cj2 



dC 

d(AM) 



m J ^ 



zitt. 



de 



70^ 



(100) 

(101) 

(102) 



for the uncertainty in local position measurement. Similarly, it becomes evident that 



V{T)\' 



mm 

leading to 



Vaidya 



2 

Vaidya 



'PiB){T) 



18 a 



2 ,.,2 



zitt. 



d(AM) 



de 



2v^ 



(AP) 



Vaidya 



Vaidya 



1 



d(AM) 



dC 



3 a 



Ar2 



r2 cj2 



d(AM) 



de 



(103) 

(104) 

(105) 



for the corresponding uncertainty for relative momentum. 

Therefore, the combination of (102) and (105) leads to the gravitationally- modified 
relative position-momentum uncertainty relation 



(AX) (AP)^^,,^^ 



1 - 



a 



zitt. 



d(AM) 



de 



256 
T57r 



+ 3i/3o)7o' (AX(B)) (AP(B)) , 

(106) 



whose diagonalized form, in terms of the chiral representation for the gamma matrices 
[16], is 

d(AM) 



(AX) (AP) 



Vaidya 



diag. 



a 



/256\2 



zitt. 



dC 



(AX(s)) (AP(^)) , 

(107) 
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where the upper sign in (107) refers to the right-handed spinor, while the lower sign refers 
to the left-handed spinor. It is interesting to note that when u^itt, 0, the absolute 
magnitude of (AX) {^'P)yaidyg, becomes infinitely large, while cuznt. oo reduces (107) 

to (AX(b)) (AP(b)). 

Assuming the upper sign in (107), the choice of « = +1 shows that infalling 
radiation serves to reduce the uncertainty relation computed in fiat space-time, though 
|d(AM) /d^l <^ 1 for realistic spherical mass accretion rates involving astrophysical 
black holes, due to the Eddington luminosity hmit [32]. As such, the gravitational 
contribution to the position-momentum uncertainty in measurement is negligibly small 
for this circumstance. However, it is unclear if this condition still applies where it 
concerns microscopic black holes, especially if r reduces to the Compton wavelength 
scale or smaller. 

Finally, to ensure that the Heisenberg uncertainty relation is satisfied, the prefactor 
in front of (^AX(b)) (^A'P(s)^ in (107) must remain nonzero. Given (97), this condition 
implies that 

d(AM) 



r > 3Mo 



1± W9 + 



/256V " 



(108) 



Again, assuming the upper sign and infalling radiation for (108), it is interesting to 
observe a growth in the lower bound for r due to gravitational and Zitterbewegung 
effects. Since r = SMq is the (unstable) photon orbit in Schwarzschild space-time, it 
is not surprising that the lower bound will grow due to mass accretion from infalling 
radiation. However, its explicit dependence on the Zitterbewegung frequency is unusual, 
a byproduct of the quantum violation of the weak equivalence principle. 

5.2. Time Evolution of the Canonical Momentum due to Zitterbewegung 
When (78) and (79) are applied to the canonical momentum, the result is 

-^^^ -Q(r)/pf(r), (109) 

= - Q(T)/P«(T) + V,Ho{T) , (110) 



dP^-\T) 



where 



0(r)/ ^ [(vanAW) 5/ - (V. lnA^'^)) 5^] cx\,){T) . (Ill) 

There are two very important points to note concerning (111). First, the fact that 
Q{T).j'' 7^ in curvilinear co-ordinates indicates that the standard free-particle approach 
only applies for Cartesian co-ordinates, which implicitly requires strictly inertial or 
rectilinear motion for the spin-1/2 particle. For non-inertial motion with rotation, a 
more complicated time evolution applies. The second point is that (111) is periodic due 
to Zitterbewegung contributions from a(Q-)(T). As a consequence, it follows that 

f7(r + 27rMitt.)/ = f^(T)/, (112) 
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which satisfies the conditions for Floquet's theorem to apply [33, 34]. This leads to the 
free-particle solution of the form 



where 



nf\T + 27i/u, 



uf\T) 



With explicit initial conditions, the free-particle solution for pf\T) is 



Pf (T) ^ Pf (0) + e^'^^»*-^nf (T) - nf (0) , 



>(0) 



(0), 



r(o), 



(113) 

(114) 

(115) 



where Re(A) > in (115) indicates an instability in the solution space for certain values 
of T and/or regions of the curved space-time background where this condition may be 
satisfied. 

To solve for P^^\t) from (110), it is possible to multiply from the left by a time- 
dependent invertible matrix ii{T)j'', such that 



/x(T)/'(V,ffG(T)) = — [^(T)/P«(T) 



(116) 



(117) 



Then so long as a solution exists for the matrix differential equation 
^^-/.(T)/Q(T),'= = 0, 

which is possible to obtain using Floquct theory, it becomes evident that lJ^{T)j^ is an 
integrating factor for (116), yielding the general solution 



/ fiiTy {V,Hg{T')) dT' + i,{0)i' Pi^)(0) 
Jo 



(118) 



for the first-order perturbation due to the interaction Hamiltonian Hq. 



5.3. Time Evolution of Covariant Spin due to Zitterbewegung 

The final set of computations concern the time evolution of the Pauli-Lubanski vector, 
the operator describing covariant spin, involving Zitterbewegung in a curved space-time 
background. For spin-1/2 particles, the Pauli-Lubanski vector in a local Lorentz frame 
is [11, 12] 

W{T) = - \ e^^fi^ a'^f'iT) P^{T) , (119) 



where a°'^{T) can be expressed as 

(T°^(r) = 2i d^^'o cc^{T) - (T^{T) . 
Applying (78) and (79) to (119) eventually leads to 



(120) 



dT 



m 



s%,ja'^,)(T)Pl,^{T) /3(o)(T) 
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+ 



h 



dTyfi)(T) 
dT 



1 



(121) 



- e%j^ {cos (u;,itt T) [^(T) a^o - 



+ sin(o7,itt.T) X»^(T)}P^o)(T)/?o 
4 



+ 
+ 
+ 



^zitt. 

i 

<^zitt. 



[cos (u;,itt.r) + sin (u;.itt.r)] ^'^f, r^,]^ C'=(r) Pj,^ (T) 
J'(T) sin (cj^itt.T) £:^oj7"o-P(o)(^) 



— rjj, - S^j S%) e'\i e'{T) PJ,^ (T) {i + Po) (122) 

l^zitt. ^ ^ 

for tiic zeroth-order and first-order time derivatives of W^{T) pa W^q^IT) + W'^^^{T), 
where 



2^1 



P^k\T),PP{T) 

= s"^'^' (Vfe In A(')) Pf\T) (123) 

in (121) is the non-inertial dipole operator [11, 12, 13], which identically vanishes for 
momentum operators in Cartesian co-ordinates. Because of the explicit and widespread 
dependence of P^q)(T') in both (121) and (122), it is possible to foresee instabilities 
emerging in the time evolution of W'^{T) due to instabilities in (115). 



6. Conclusion 



This paper explores the consequences of spin-1/2 particle Zitterbewegung in a general 
curved space-time background using Fermi normal co-ordinates, where the spatial 
vectors orthogonal to the classical worldhne are treated as quantum fluctuations 
X|^^(T) in a local frame, according to the formalism first introduced by Barut and 
Bracken [14]. An overarching conclusion is that the associated Dirac Hamiltonian Hq 
reveals a quantum violation of the weak equivalence principle due to Zitterbewegung, 
confined to a worldtube of radius u;^^\ . Apart from one anti-Hermitian term in (75) 
that is coupled to ^ R^^jk,o{T), all the gravitational contributions in Hq are dependent 
on the Ricci curvature tensor and its time derivative in the Fermi frame. In the classical 
hmit as ^ — > 0, it follows that Hq — > 0, as expected. It should also be noted that, while 
the definition of quantum operators employed in this paper is sufficient to demonstrate 
the possible existence of Zitterbewegung effects in a curved space-time background, a 
more mathematically rigourous definition may be required to illustrate the full extent 
of these effects. This is a worthwhile consideration for future research on this topic. 
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When applied to the special case of quasi- circular orbital motion in the presence 
of a spherical black hole in radiation, as described by the Vaidya metric, it is shown 
that uncertainties in the simultaneous measurement of the particle's locally determined 
position and momentum are modified due to the coupling of Wznt. with the black 
hole's mass rate-of-change |d (AM) /d^\, with interesting physical consequences. By 
further allowing a description of the problem in terms of general curvilinear co-ordinates 
to accommodate non-inertial motion along the worldline, the purely fiat space-time 
treatment of this problem becomes generalized by showing that the canonical momentum 
is subject to the Floquet theory of differential equations involving oscillation frequency 
'A^zitt., with potential stability implications due to Zitterbewegung. 

There are numerous potentially significant consequences and future directions worth 
exploring from the main expressions presented in this paper. For example, the fact 
that Hq is explicitly time-dependent due to Zitterbewegung indicates the possibility of 
identifying transition rates for local energy states about the spin-1/2 particle's classical 
worldline [34] . While the gravitational and Zitterbewegung contributions do not survive 
for a c/asszca/ vacuum space-time background, where ^R^i,{T) = identically, this may 
not necessarily be true if zero-point fluctuations and vacuum polarization effects [35] due 
to the Casimir effect, for example, appear as a consequence of virtual pair-production. 
Given that QED predicts the existence of vacuum bubbles that have no observable 
effect in purely fiat space-time, it is reasonable to consider whether this property holds 
true in the local presence of space-time curvature. Finally, the results presented here 
may have potential cosmological implications in terms of the very small value for the 
(observed) cosmological constant, with the perceived global acceleration of the Universe 
as an observable effect. These and other possibihties may be explored at a later date. 
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